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Abstract 

In this note, we give a probabilistic interpretation of the Central Limit 
Theorem used for approximating Gaussian filters in 

It was shown in ^ how a certain "radial" form of the Central Limit Theorem 
(CLT) could be used to approximate isotropic Gaussians on the plane. The main 
idea was to approximate the Gaussian using box distributions that had been 
uniformly distributed over the half-circle. The term "radial" was used here to 
highlight the fact that the Gaussian filtering was achieved by convolving the 
image using box distributions along radial directions (and parallel to it) . 

The idea was extended in two different directions in fll . In one direction, it 
was argued that, by adjusting the widths of the box function along each radial 
direction, one could approximate anisotropic Gaussians. In particular, for the 
special case of the so-called four-directional box splines, a simple algorithm was 
developed that allowed one to control the covariance by simply adjusting the 
widths of the box distributions. In a different direction, an algorithm for space- 
variant filtering was developed which, unlike convolution filtering, allowed one 
to change the shape and size of the box spline at each point in the image. This 
was also done at the expense of just 0(1) operations (independent of the shape 
and size of the Gaussian) using a single global pre-integration followed by local 
finite-differences; cf Algorithm 1 in 

A probabilistic interpretation of this result is as follows. Let X be random vector 
on the plane that is distributed on a line passing through the origin (e.g., one of 
the coordinate axes). Thus, X is completely specified by a probability measure 
li{t) on the real line. Suppose that 



Guided by the idea behind Theorem 2.2 in [[U, we can then conclude the 
following. 




For < < TT, let us denote the rotation matrix on the plane by 
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nRotations = 10; 




nDraws = 10000; 




rotation = pi*(0: 


1: (nRotations — 1 )) /nRotations ; 


Zl = zeros (nDraws, 


1) ; 


Z2 = zeros (nDraws, 


1) ; 


for n = 1 : nDraws 


% outer loop (realizations) 


XI = 0; 




X2 = 0; 




r = sqrt ( 12 )* (rand (nRotations, 1 ) — . 5*ones (nRotations , 1 )) ; 


for k = 1: nRotations % inner loop (rotations) 


XI = XI + 


r (k, 1) *cos (rotation (1, k) ) ; 


X2 = X2 + 


r(k,l)*sin{rotation(l,k)); 


end; 




Zl(n,l) = sqrt (1/nRotations) *X1; 


Z2(n,l) = sqrt (1/nRotations) *X2; 


end; 





Figure 2: The MATLAB simulation of the Central Limit Theorem. 



Theorem 1 (Radial CLT). For any integer N, fix Ok ^ {k - 1)tt /N for I < k < N. 
Let Xi, X2, . . . , Xjv be independent and identically distributed copies of X, and 
set ^ 

Zn = —^ {TZg^Xi + TZe2X2 + ■ ■ ■ +TZ0j^Xn) . 

Then the random vector converges in distribution to the standard normal 
distribution as N gets large. More precisely, for any Borel set B on the plane, 

lim ProbfZAT e S) = — / exp f-^^^ dx. 

Indeed, the proof of Theorem 2.2 in [fTT] shows us that that the characteristic 
form of Zjv converges to that of the normal distribution as N gets large. The 
weak convergence of the distribution of thus follows from Levy's theorem. 

We performed a simple MATLAB simulation for the case where X is \/l2 times the 
uniform distribution (zero mean, unit variance) . The results of the simulation 
when iV = 10 is provided in Figure [2] 
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